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ABSTRACT 



A Markov type model for studying rank distribution in 
an hierarchical organization is examined. Various sets of 
rank distributions are defined and their properties discussed. 
A computer aided test is developed for testing a given rank 
distribution to determine if it is an element of a particular 
set. For a distribution not in the Steady State set, a test 
is developed to determine if the distribution can be returned 
to in m transitions or steps. 
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I. INTRODUCTION 



Several authors have studied models for analyzing rank 
distributions in an hierarchical organization. Bartholomew 
[1967] uses Markov models to a great extent in his analysis 
of various social systems. Thonstad [1968] developed a 
Markov type model for use with the Norwegian educational 
system. More recently, Branchflower [1970] utilized a Markov 
type model to examine the promotion and appointment policies 
of a university and the resulting rank distributions of 
faculty. Branchflower assumed the promotion policies were 
fixed and used his model to predict the faculty distribution 
in the next time period using the current faculty distributions 
and actual appointments to each rank. In his analysis, 
Branchflower defines three sets of rank distributions: (1) an 

Attainable set, (2) the Maintainable set, and (3) the 
Containment set. An Attainable set is a set of rank distri- 
butions which can be reached from some previously attainable 
set. The Maintainable set is better known as the Steady 
State set of rank distributions. The Containment set is the 
set of rank distributions which. can be returned to infinitely 
often by means of a sequence of appointment policies. 

The analysis in Branchflower [1970] leads us to pose 
several questions about the rank distributions given some 
promotion policy: 

1. Can the present rank distribution be maintained for 
all time by means of appointment policies alone? 
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2. Can another given rank distribution be attained from 
the present distribution in a finite time period, or can the 
desired distribution never be reached? 

3. What are the rank distributions which can be returned 
to at some time? 

In this paper, our purpose is first, to define explicitly 
the sets of rank distributions of interest associated with 
the model, and to demonstrate various properties of these 
sets. Second, we shall use these definitions and properties 
to develop test procedures to determine if a given distribu- 
tion is an element of a particular set. This analysis will 
provide answers to questions 1 and 2, and will permit us to 
give crude bounds as an approximation to answer question 3. 
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II. THE BASIC MODEL AND NOTATION 



Consider a discrete time Markov model with n+1 states 
where state j is the rank of a member of an hierarchical 
organization. The set of states can be partitioned into two 
subsets: the first n states are transient states and are 

called active ranks; the second subset is the n+l s *- state 
which is an absorbing state. Let p^j be the fraction in 
state j at time t+1 that were in state i at time t. Then 



n+1 

p • • 1 , 1 — 1 , 2 / . . . , n+ 1 . 

j=l D 



We represent these fractions in a transition matrix 



Q 



P 

0 



w 

1 



where P is an n x n matrix, [p^j], i , j=l , 2 , . . . ,n , and w is 
an n dimensional column vector called the loss vector with 
Wj_ = p^ n+ i* P is called the Promotion matrix and describes 
the transition process of the active ranks in the organization. 

n 



It has the property that 



j=l 



ij 



< 1, i=l , 2 



n with strict 



inequality holding for at least one row. 

We assume that the organization is in an equilibrium 
condition, and thus keep the total number in the organization 
constant. We normalize by assuming this number is one. 

Let x(t) be an n dimensional row vector, x(t) = (x-^(t), 

i 

X 2 (t) , . . . , x R (t) ) where x^ (t) is the fraction in the ith rank 
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at time t . 



Note that E x^ (t) = 1, (t) 0. Let f (t) be 

i=l 

an n dimensional row vector, f (t) = (f ^ (t) , . . . , f (t) ) where 
fi (t) is the fractional number of new appointments to rank 
i at time t, f^(t) >_ 0. 

The rank distribution at time t+1 is dependent on the 
rank distribution at time t, the Promotion matrix, and the 
appointment policy at time t+1. In vector notation, 

x ( t+1 ) = x (t ) * P + f(t+l). (1) 

The requirement that the total number in the organization be 
constant requires that the total fraction appointed to the 
active ranks be equal to the total fraction lost to the 

" v o' n 

absorbing state: Z f . (t+1) = x(t)»w. (2) 

i=l 

We use the following notation throughout the remainder 
of the paper: e is an n dimensional vector with all components 

equal to 1; ej is an n dimensional unit vector with the 
component equal to 1. It is apparent by the context whether 
they are row or column vectors . 



III. SET DEFINITIONS 



The set of all rank distributions, F, is defined as 
F = {x; x*e = 1, x ^ 0} 

i 

Note that the time notation is omitted when it does not affect 
the definition. 

An Attainable set is a set of rank distributions which 
can be reached in one step from a previous set of rank distri- 
butions. Let m be the number of steps (time periods) from 
the set of all rank distributions, F, then 

A^j = {x; x >_ z*P, x e F, z e A m -i^ > m = 1 ,2,... and 
A q = F. 

Consider the sequence of Attainable sets, {A^}. If the 
sequence converges to a set C we call the limiting set the 
Containment set. Mathematically, we have 

C = lim Ajjj 
m -»■ 00 

= lim {x; x^z*P, x £ F, Z £ A m _^} 
m -»■ 00 

= (x; x >_ z*P, x e F, z e C}, if the limit exists. 

This definition is somewhat unusual in that the Containment 
set is defined in terms of itself. 

The set of all distributions which can be maintained in 
successive time periods is called the Steady State set, M. 

Let f (t) = f for every t, then since x(t+l) = x(t), from 
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(1) we have x(t) = x(t)*P + f 
x (t) • (I-P) = f > 0 



so we define 

M = {x; x* (I-P) 0 7 x e F} 

If a rank distribution can be reached from the Steady 
State set, M, then we call a set of such distributions a 
Reachable set. All distributions which are reachable from M 
in m steps are in the set B m , where 

B m = {x; x j> z*P, x e F, z e B m _i) / m = 1,2,... and 
B q = M. 

Now consider the sequence of Reachable sets, {B^}. If 
this sequence converges to a set B, we call the limiting set 
the Reached set. Mathematically, we have 

B = lim B m 
m -*■ 00 

= lim {x; x>^z*P, x e F, z e B rn-1^ 
m ■+ 00 

= {x; x z‘P, x e F, z e B) , if the limit exists. 

Thus by their definitions, we see that the Containment set and 
the Reached set must be the same set when they both exist. In 
section IV we prove they always exist. 
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IV. SET PROPERTIES 



In this section, we list several properties of the sets 
defined in section III without proofs. We then prove several 
lemmas and two Theorems which establish the existence of the 
Containment set and Reached set. 

PROPERTY 1: The set of all rank distributions, F, is 

convex and closed. 

PROPERTY 2: Any Attainable set, A Jn , is convex and closed. 

PROPERTY 3: The Steady State set, M, is convex and 

closed. 

PROPERTY 4: Any Reachable set, B , is convex and closed. 

J m 

LEMMA 1: An Attainable set, A^, is a subset of every 

previously attainable set. 



Proof : 


F = Aq = {x; x*e 


= 1, x > 


0} 






A^ = {x; x _> 


z *P, x e 


F, z 


£ Aq} 


Pick x 


£ A^ ^ x e Aq 


there fore 


A 1 £ 


> 

o 

II 


Assume 


Am-l £ %-2 5 • 


Ul 

1 — 1 
< 

Ul 

• 

• 


Aq. 




Define 


A m (z) E {x; x >_ 


z*P, x e 


F} , 


z e A m _! so 



A m " U V z) ' pick z * e *m-l» then since \n-l = *m-2' 

z e Vl 

A m (z*) = {x; x _> z**P, x e F, z* e A Jtl _ 2 } / and 
A m <z * ) = ^m-l * 

So % ‘ U *nif z >- c Vl- 

z e A , 

m-1 

LEMMA 2: The Steady State set, M, is a subset of any 

Attainable set, A . 

m 
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Proof: M = {x; x _> x*P, x e F}. M £ = F by definition. 

Assume M c A , . Pick z * e M ^ z * e A and z* > z**P. 

m~± m—± — 

Therefore, z* e A m by definition of A m , and M £ A m . 



LEMMA 3: Any Reachable set, B m , contains every previously 

reachable set. 

Proof: M = Bq = {x; x _> z*P, x e F} . Trivially, Bq £ B^. 

Assume that B m _2 c B m-1 w here B m -i = {x; x >_ z*P, 
x e F, z e B m _ 2 ). 

Define the following sets: B m _-^(z) = (x; x ^ z*P, 

x e F}, z e B m _ 2 , an< ^ B m (z) = {x; x z*P, x e F}* 1 , 

z e Bm-i. Pick z* e B m _ 2 . Since B m _ 2 £ B m _ 1 , z* e B m _ lf 

so that B m _^(z*) = B m (z*). Since z* is arbitrary, and 

since B m _ 2 £ Bm-i, B^ = U B m _ 1 (z)= U B m (z) 

z e B m _ 1 z e B m _ x 

= B . 
m 

LEMMA 4: If w^ > 0 for at least: one i and n > 1, then 

the Steady State set, M, contains an infinite number of points. 
Proof : 



(i) Non-empty . Assume that M is empty. Then there 
does not exist an x such that x* (I-P) >. 0* Because 
n 



I 

j=l 




< 1 / 



i=l, 2, 



,n with strict inequality holding 



for at least one i, (I-P) - -*- exists with all elements 
non-negative. Let N = (I-P) -1 . If M is empty, then 
there does not exist an x such that x = f»N ■> 0. But 
n 

Xj = E f^*n^j and since f^ 0 and n^j >_ 0 for every 
i,j, Xj > 0 and therefore M is non-empty. 
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(ii) Infinite Interior. Assume that M contains only 



one point, then x(t + 1) = x(t) for every t and any 
appointment vector. We rewrite equation (1) as x(t) = 
x(t)*P + f(t + 1) = f (t + 1) • (I-P) “-*• . We also require 

that f (t + 1) *e = x(t)*w. We can construct the 

appointment vector in the following manner: let 

f(t + 1) = cue. + a^e. +...+ a e where Za = x(t)*w. 

1 1 22 nn ii 

Then f (t + 1) = (a^ , , . . . / a n ) and 

x (t) = (a 1 ,a 2 , . . . ,a n ) • (I-P) (3) 

Consider the appointment vector f (t + 1) when 
cc^ = x(t)*w and a 2 = a 3 = = a n = ® ' then x (t) equals 

x(t)*w multiplying the first row of (I-P) - '*'. Since 
any f which solves (3) generates a point in M, consider 
the appointment vector f (t + 1) when 012 = x(t)*w and 
a l = a 3 = ••• = a n = 0, then x(t) equals x(t)*w multiply 
ing the second row of (I-P) - -*-. Since x(t) is one point, 
and x(t)*w is a scalar, the first two rows of (I-P) - -*- 
are identical. This implies that r((I-P) - -*-) < n and 
that ( (I-P) ”'*') "I does not exist which is a contradiction 
Thus since M is non-empty, M contains more than one 
point. Since M is convex and closed, it contains an 
infinite number of points. 

THEOREM 1: The Containment Set, C, and the Reached Set, 

B, exist, are identical, and non-empty. 
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Proof: Consider the sequence of sets {A^} . By Lemmas 

1 and 2 , M £ . . . c E A s . . . £ A n c F . By Lemma 

m m-1 1 

4 M ={= <j) => A m : f<j) f n = l,2 / ... . Thus we have a mono- 
tone decreasing sequence of sets which is bounded below, 
therefore the limiting set does exist and is non-empty. 

Now consider the sequence of sets (B m ) where Bq = M. 

By Lemma 3 , M £ B ft £ B, £ ... SB , £ B £ . . . £ F . 

u u ± m-l m 

Thus we have a monotone increasing sequence of sets 
bounded above by F. Therefore a limiting set exists. 

As noted in section III the Containment set, C, and 
the Reached set, B, are the same set. 

THEOREM 2: The Steady State set, M, any Attainable set, 

A m , the Containment set, C, any Reachable set, B m , and the 
set of all rank distributions, F, are identical if and only 
if the Promotion matrix, P, is diagonal. 

Proof: Since M £ B m £ C £ Aj^ c F, we need only show that 

M = F. 

(i) Assume that P is diagonal. We can write the 

distribution vector as x(t) =ae,+ae +...+ a e , 

1 1 2 2 n n 

where a. >0. If Ja. =1 then x(t) e F. Rewrite 
i — ii 

equation (1) x(t + 1) = x(t)*P + f(t + 1). We require 
f (t + 1) *e = x (t) • (I-P)*e . Since P is diagonal, x(t)* 

(I-P) is a suitable choice for f (t + 1), (f^ (t + 1) ^0). 

So we can now write x(t + 1) = x(t)*P + x(t)*(I-P) = x(t). 
Now x(t + 1) = x(t) and x(t) e F =* x(t) e M, F £ M. 

And since M £ F, M = F. 
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(ii) Assume M = F. x(t + 1) = x(t) for all t. 

Assume that P is not diagonal. x(t) = x(t)*P + f (t + 1) . 
since x(t) e F, we can write x(t) = a^e^ + & 2 e 2 + ••• + 
a n e n' «i 1 0, =1. In particular, let = 1, 



l 



= 0, i 1. So x(t) 



( 1,0 



0) 



Pll P 12 



= a^e^ and (1,0,..., 0) = 

‘ * * Pin" 



+ f(t + 1) 



_ p nl p n2 “• P nn 

- 7 P 12 7 • • • > P}.n^ ^1 7 -^2 7 • • • 7 * This requires 

that the following equations be satisfied: 



f 2 + Pl2 = 0 

• • • 

• • • 

• • • 

f n p ln “ ® 

Since f^ > 0, p.. > 0, we have f . = p. . = 0, j I 1. We 
repeat this procedure for otj = 1, = 0 , k=fj, j = 

2,3,...,n. Thus p^j = 0, i =}= j and P is diagonal. 
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V. EXTREME POINTS OF THE VARIOUS SETS 



The characterizations of the various sets in the model 
in the above section do not describe them in a useful manner. 

If we knew the extreme points of all of the sets, then we 
would have almost all of the information about the sets, 
since any interior point of a set is simply the convex combin- 
ation of its extreme points. 

For two sets in particular, the extreme points are easily 

obtained. The set of all rank distributions, F, is the 

fundamental simplex in E n with extreme points, ej_, i = 1,2, ...n. 

The subset of F which is most easily described is the Steady 

State set, M. Each vector in the set satisfies the inequality 

x x*P. This can be written in equation form as x* (I-P) = f 

where f is the vector of appointments to the different ranks. 

We define N = (I-P) - -*- so we write x = f*N. Since we require 

x*e = 1 and since f can be written f = Za-^ (x* (I-P) *e) , ^ 0, 

i 

= 1, we write the equation in an equivalent form as x = f*B 
i 

N- • 

where f = Za-e. and B- . = . That is, the extreme points of 

i 1 1 13 Pij 

3 

the Steady State set are given by the rows of (I-P) - ^ after the 
rows are normalized. 

The problem of determining the extreme points of the other 

sets in the model is much more difficult. If we are working in 

three dimensions, it is not too difficult to calculate all of 

the extreme points of one of the sets, A or B , for a small m. 
c mm 

• I 



20 



Each set in the model has the property that each point in the 
set generates a subsimplex in the fundamental simplex. The 
extreme points of the generated set are found from the convex 
hull of the subsimplices generated by the extreme points of 
the original set. In three dimensions we use a geometrical 
method to determine the extreme points of the convex hull of 
the generated subsimplices. If we know the extreme points of 
for some m, then we project A m onto a plane generated by 
two unit vectors. Then we calculate the extreme points of 
the subsimplex generated by each extreme point of A m and project 
these onto the plane. Then we simply test linear combinations 
of the generated extreme points to determine that there are 
no points which lie beyond the linear combination in each of 
the unit directions. The points which satisfy this criterion 
are the extreme points of the set A m +]_. A similar procedure is 
used for B m . 

When P is upper triangular, some of the extreme points of 
A m are also extreme points of A m+ ^ (e.g. e n is an extreme 
point of Aj^ for every m.) These remain extreme points for all 
succeeding sets generated from A m+ ^. Similarly, B m has 
extreme points common to all Bj^ which are extreme points of 
A m as well, and therefore are extreme points of the Containment 
set, C. Figure 1 depicts several attainable sets drawn on 
barycentric coordinates for a three dimensional rank 
distribution where 



21 




22 



.5 



.4 



0 



P = 



0 

0 



.6 

0 



.3 

.7 



Table I contains the coordinates of the numbered points on 
Figure 1. The heavy line in the figure represents the known 
portion of the boundary of the Containment set after calculating 
A-]. The portion of the boundary between points 9 and 10 
changed very little between Ag and A^ so the boundary of A 7 
is probably a good approximation of the boundary of the 
Containment set in this case. 



TABLE Is POINTS ASSOCIATED WITH FIGURE 1. 



Point 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Rank 


1 


0 


.3 


.39 


.393 


.3723 


.35289 


.18615 


.176445 


.0882225 


0 


1/3 


2 


0 


0 


.12 


.228 


.2940 


.32532 


.49206 


.500706 


.5355714 


.5 


1/3 


3 


1 


.7 


.49 


.379 


.3337 


.32179 


.32179 


.322849 


.3762061 


.5 


1/3 



It is not clear how this graphical projection method can 
be extended to higher dimensions to obtain the extreme points. 
Nor is it clear, based on the above example, whether or not 
the Containment set has a finite number of extreme points. 
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VI. THE DEVELOPMENT OF SEVERAL TESTS 



Since it is not feasible to analytically determine the 
extreme points of all of the sets, we must devise some method 
to answer our original questions. The answers to those 
questions depend on knowing the extreme points of all of the 
sets. To provide partial answers to the questions, we develop 
test procedures to determine if a given vector is an element 
of a particular set, and to determine if a given vector can 
be returned to in m steps or reached from another specified 
vector in m steps. We know that a given vector cannot be 
returned to at any time if there exists an i such that the 
test vector is not in Ap 

Each transition from one set to another requires that a 
certain set of inequalities be satisfied. When we consider 
the system of inequalities that account for all of the transi- 
tions, we need only test for feasibility. If no feasible 
solution exists, then the test vector fails the test. Several 
specific tests are developed in the next sections. 

A. A TEST FOR 

m 

In this section, we develop a test to determine if a given 
rank distribution is an element of the set B m . Recall that 

B m = {z m ; z m 1 z m-l* P ' z m e F ' z m-l e B m-1 } for m = 1/2,... 
and that B m is the set of rank distributions which is reachable 

in m steps from the Steady State set. Row vectors have been 
used in all developments. To simplify notation, the test 
inequalities will be specified in tableau form and vectors 
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will be column vectors. Transposition is indicated for the 
promotion matrix where it appears. 

We now write the sequence of inequalities which must be 
satisfied if the vector X is in the set B m as 



z 0 Z 1 z 2 
(P T -I) 



e 



e 



z 0 z 
m-2 m-1 



< 0 

1 0 

< 0 




< 0 



< X J 



( 4 ) 



= 1 



e 



= 1 



= 1 



e = 1 



The constraint matrix appears deceptively smaller than it 
is. Since each is an n vector, the dimensions of the matrix 
are mn + m + n rows and mn columns . Even though the constraint 
matrix is very sparse, we were unable to apply any algorithm 
which would process the problem more efficiently than the 
Simplex algorithm. 
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It is clear that a feasible solution to (4) is simply a 

l 

Pha^e I Simplex solution. The method which we use to obtain 
this solution is to add an artificial variable to each of the 
feasibility constraints, the last m rows. We then minimize 
the sum of the artificial variables. A feasible solution to 
the new system, (4) plus the artificial variables, always 
exists. The result of the test is determined by the values of 
the artificial variables. If all of the artificial variables 
are equal to zero, then (4) has a feasible solution and we 
conclude that the test vector, X, is an element of B m . If 
any of the artificial variables are not zero, then a feasible 
solution to (4) does not exist and we conclude that X is not 
in B m . 

B. A TEST FOR A m 

In this section, we develop a test to determine if a 
given rank distribution is an element of A m , i.e. if it can 
be reached in m steps from some point in F. Recall that 

% = {z m'* z m 1 z m-l * P ' z m e F ' z m-l e ^n-l* for m = 1/2,... 

We write the set of inequalities as 
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z 0 Z 1 z 2 * 2 in- 2 z m-l 

P T -I < 0 

P T -I < 0 



e 



e 



pT 



e 



-I £ 0 

p T < X T (5) 



e 



= 1 

= 1 
= 1 

= 1 



The constraint matrix has mn + m rows and mn columns. We 
now test (5) for a feasible solution in the same manner as 
the test for B m . The results of the test are interpreted in 
the same manner as those for B^. Note that the constraint 
matrix in (5) is the same as that in (4) without the first n 
rows . 



C. A TEST FOR RETURN IN m STEPS 

In this section, we develop a test to determine whether 
or not a given rank distribution, X, can be returned to in m 
steps. If X is an element of the Steady State set, we know 
that it can be returned to in one step, i.e. it can be 
maintained forever. So we are interested in a test for m 
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greater than one. Since there are m transitions, there must 
be m-1 intermediate vectors because the initial and terminal 
vectors are both X. The first intermediate vector, z^, must 
satisfy z^ >L X‘P. The second intermediate vector, z 2 , must 
satisfy Z 2 _> z^*P. This process continues until the last 
intermediate vector, z m _i/ must satisfy z m _-^*P _< X. We write 
the sequence of inequalities in tableau form as 



z 



1 



Z - Z -a 

m-2 m-1 



I 

pT 



-I 

„T 



P T. X T 



< 0 



< 0 



,T 



-I <0 



PT 



1 X 

= 1 

= 1 

= 1 
= 1 



T 



( 6 ) 



Again we notice that there is a great similarity between 
the form of this tableau and those for the tests for A jn and 
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B m . The first n inequalities in (6) are different from those 
in (4) , as well as there being one less intermediate vector. 

The constraint matrix in (6) has mn + m + n rows and mn - n 
columns . 

As with the tests for A m and B^, any method of finding 
a feasible solution to (6) is acceptable. One method would 
be to add artificial variables not only to the last m rows 
of (6) , but also to the first n rows, then minimizing the 
sum of the artificial variables. If all of the artificial 
variables are zero, then we conclude that the rank distri- 
bution, X, can be returned to in m steps. The method used 
in section VII involved adding artificial variables only to 
the last m rows of (6) , then minimizing the sum of the artifi- 
cial variables. This system may not always have a feasible 
solution. The test requires that if the system has a feasible 
solution, then if all of the artificial variables are zero, 
we conclude that X can be returned to in m steps. 

D. A TEST FOR CHANGING DISTRIBUTIONS IN m STEPS 

We now develop a test to determine, given some distribu- 
tion vector, X^, whether or not we can reach another specified 
distribution vector, X£/ in m steps. This is nearly identical 
to the problem discussed in the above section, i.e. we have a 
fixed initial vector and a fixed terminal vector. The con- 
straint matrix is the same as (6) but now the right hand side is 
[P T -XT, 0,..., 0, X*, 1,..., 1]. The test procedure is then 
the same as the test for return in m steps. 
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E. A TEST FOR EXTREME POINTS OF THE CONTAINMENT SET 

The similarity between the tests for A^ and B m noted in 
sections A and B above permits us to develop a procedure with 
which we can determine extreme points of the Containment set, 

C. 

T 

Let Gji = [ (I-P ), 0,..., 0] be an n by mn matrix and let 




e 

-e 



e 

-e 



e 

-e 

be an mn + 2m by mn matrix. Now let = [0, . . . ,0,X T / 1, . . . ,1] 

T rn 

be an mn + 2m by 1 column vector and Z = [z , ... ,z x ] be an 

0 m-1 

mn by 1 column vector. Then the test for A^ reduces to finding 
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a feasible solution to 




( 7 ) 



If there is a Z* which satisfies (7) then the test for 
would seem to result in testing to see if G^*Z* 0. However, 

there is no guarantee that Z* is unique and there could 
exist a Z' which solves (7) and also solves G^*Z' >_ 0 for the 
same test vector X. So this would be a sufficient test for 
inclusion in B m , but not necessary. The only necessary and 
sufficient test is that described in section III. A. 

Now let H-^ = [ (I-P T ) , 0 , . . . , 0] be an n by mn + n matrix 
and let 



P T -i 



P 



T 



-I 



I 




e 

e 



H 



2 



e 

e 



e 

e 



e 

e 
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[ Z /->/•••/ 2 ] be 
u m 



be an mn + n + 2m by mn + n matrix. Let Z' = 

an in + n by 1 column vector and let = [0,...,0,1,...,1] 

be an mn + n + 2m by 1 column vector. This amounts to letting 

the test vector, xT become an unknown, z T . Then all of the 

m 

basic feasible solutions to I^’Z' _< X^ are extreme points of 
A^. If a point is an extreme point of A m and B m simultaneously, 
then it is an extreme point of the Containment set, C. Thus 
we need only test each basic feasible solution to see if 
H^*Z' ^0. If it is, then Z' is an extreme point of C. 

Due to the computational difficulties associated with 
this test, it is presented for theoretical interest only. 
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VII. A CASE STUDY 



The tests developed in section VI above were applied to 
faculty distributions in the College of Engineering at the 
University of California at Berkeley based on data from 1960 
to 1968. The bulk of the data is in Branchflower [1970] and 
only that which is relevant is reproduced here. 

We used the Mathematical Programming System on the IBM 
360/67 computer at the W. R. CHURCH Computer Center, Naval 
Postgraduate School to conduct most of the tests. This is 
really a trivial exercise once the data cards are prepared 
for the control programs. One FORTRAN IV program was written 
to generate the data cards for the test for B m . The same 
cards were used to test for Ajjj after removing the cards per- 
taining to the first n rows of the constraint matrix. Another 
FORTRAN IV program was written which was used to generate 
data cards for testing return in m steps. Again, hand modifi- 
cation of the right hand side permitted the same cards to be 
used to test for changing distributions in m steps. 

One other FORTRAN IV program was written which was used 
to determine the extreme points of the steady state set. It 
was also used to test to determine if a given vector was an 
element of the Steady State set. 

All of the tests were verified for accuracy by using a 
test case, namely the three state promotion matrix in section 
V and points which were in known sets from Figure 1. 
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k 



















STATE 



I 


1 2 


3 4 


5 


6 


1 


.28 .56 


.04 .04 






2 


.41 


.47 .07 






3 




.44 .37 


.14 




4 




.38 


.52 


.04 


5 






.47 


.43 



10 11 12 13 



6 

7 

8 
9 

10 

11 

12 

13 



.57 .41 .01 
.64 .34 

.67 .30 .01 

.69 .25 .02 

.71 .24 .01 

.73 .25 .01 
.77 .20 
.97 



FIGURE 2. PROMOTION MATRIX FROM 1960-1968 DATA 



TABLE II. TEST VECTORS 



\VECT0R 










STATEN. 


X].(1965) 


X 2 (1966) 


X 3 (1967) 


X 4 (1968) 


1 


.0161427 


.0154639 


.0 


.0047619 


2 


.0802139 


.0567010 


.0628019 


.0190476 


3 


.0534759 


.0979381 


.0821256 


.0809523 


4 


.0374331 


.0257731 


.0628019 


.0761904 


5 


.0534759 


.0618556 


.0483091 


.0523809 


6 


.0481283 


.0463917 


.0628019 


.0666666 


7 


.1390374 


.0979381 


.0821256 


.0952380 


8 


.1176470 


.1494845 


.1400966 


.1047619 


9 


.1122994 


.1134020 


.1159420 


.1238095 


10 


.1390374 


.1030927 


.1062801 


.0904761 


11 


.0962566 


.1185567 


.1014492 


.1047619 


12 


.0855614 


.0670103 


.0676328 


.1047619 


13 


.0213903 


.0463917 


.0717328 


.0761904 
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For the case study, we used the thirteen active states 
of the faculty distribution and the 1960-1968 promotion 
matrix from Branchf lower [1970 ] . This is reproduced as 
Figure 2. The actual distributions of faculty for 1965-1968 
were normalized and used as the four test vectors. These 
are included in Table II. 

The test vectors were tested for being elements of A m 
and B m , m = 1,2,... ,7, and were tested to determine if they 
could be returned to in m steps for m = 2,3, ..., 8 . The 
resulting location of the vectors of rank distributions on 
the fundamental simplex is indicated in Table III. 

TABLE III. TEST RESULTS 



X 1 


X 2 


X 3 


x 4 


X 1 e 


A 2 




X 3 e A 2 


X 4 E A 2 


Xy $ 


A 3 


X 2 ^ A 1 


X 3 * A 3 


x 4 4 a 3 



Recall that if e A^ then after one step, X^ + ^ _> X^*P, 
X^ + ^ is at least an element of Aj + ^. The results in Table III 
are inconsistent with this theoretical result and require an 
explanation. There are two reasons why this could have occured. 
The first reason is that the Promotion matrix was an average 
over an eight year period and did not exactly represent the 
true Promotion matrix for each year. The second reason is that 
the faculty size did not remain constant, but in fact was 
increasing. An expanding organization can appoint personnel 
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to replace the losses and then hire the additional number 
for expansion into any state. When P is upper triangular, 
the extreme points of successive A m may tend to move away 
from the lower ranks. For a three state example, consider 
an initial rank distribution of (1, 0, 0) with 



P = 



.5 .4 0 

0 .6 .3 

0 0.7 



then after one transition, the rank distribution before new 
appointments is (.5, .4, 0) . Say the number of persons in 
the organization is constant at 100, then the distribution of 
persons in the ranks is (50, 40, 0) with a loss of 10 to 
replace. Appointments could result in the distributions 
(60, 40, 0), (50, 50, 0), or (50, 40, 10) or any linear 

combination of them. If the organization expands by 10 persons, 
then the distributions could be (70, 40, 0), (50, 60, 0), or 

(50, 40, 20) or any linear combination of them. If we 
normalize, then these expanded distributions are outside A^ 
for the constant size case. The result is that the possible 
distributions for the expanding size model can be outside of 
the normal set of distributions obtained for the constant size 
model . 

We then aggregated several ranks in the model to determine 
if the aggregated model behaved in the same way. The first 
four ranks are Assistant Professors, the next three ranks are 
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Associate Professors, and the last six ranks are Full 



Professors. The resulting Promotion matrix is 



7649 


.1752 


0 


0 


.8030 


.1841 


0 


0 


.9707 



The normalized rank distributions for 1965-1968 are 
listed in Table IV. 



TABLE IV. TEST VECTORS: AGGREGATED MODEL 



^^jector 

ramkX. 


X 1 (1965) 


X 2 (1966) 


X 3 (1967) 


X 4 (1968) 


Asst . 
Prof. 


.187166 


.195876 


.207729 


.180952 


Assoc. 

Prof. 


.240642 


.206186 


.193237 


.214286 


Full 

Prof. 


.572192 


.597938 


.599034 


.604762 



These vectors were tested to determine their location on 
the fundamental simplex. The results are shown in Table V. 



TABLE V. TEST RESULTS: AGGREGATED MODEL 



*1 


x 2 


x 3 


x 4 . 


X, t A 


X e A 


X e A 


X e A 


1 9 


2 7 


3 6 


4 9 


X 1 ♦ A 10 


X 2 i a 8 


X^ £ A 7 


X 4 * A 10 
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These results using the aggregated model are similar to 
the results of the thirteen state model. We noted above that 
as faculty size increases, the rank distribution may move to 
a superset of the set containing the original distribution. 

In fact, this happens in 1966 for an increase of seven faculty 
members (3.7%), and also in 1967 for an increase of thirteen 
faculty members. In 1968, there was an increase of three 
faculty members, but the rank distribution moved to a subset 
of the set containing the 1967 rank distribution. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



We have established several properties of several sets 
of rank distribution in an hierarchical organization which 
have not been previously examined. In addition, based on 
these properties, we developed several tests to determine if 
a specified rank distribution is an element of certain sets. 

The tests were applied to data from a case study of the 
College of Engineering faculty at the University of California, 
Berkeley. 

Several questions remain unanswered and are possible 
areas for further study. Is there some neighborhood of a 
given rank distribution which can be maintained? Is there 
an efficient algorithm which will solve the test inequalities 
other than the Simplex method? Is there a method of reducing 
the number of basic feasible solutions which must be found to 
determine the extreme points of the Containment set? Does 
the Containment set have a finite number of extreme points? 
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* * * * DATA GENERATOR - TEST FOR RETURN IN M STEPS * * * * 



* * * * * * ROBERT Lo ARMACOST - - - 1 JULY 1970 * * * * * * 






THIS PROGRAM WILL PUNCH OUT A SET OF DATA CARDS FOR A 
LINEAR PROGRAM USING THE MPS PACKAGE ON THE IBM 360/67, 
THE BASIC INPUT TO THE PROGRAM IS A SQUARE PROMOTION 
MATRIX OF SIZE N, AND K VECTORS WHICH ARE TO BE TESTED 
TO DETERMINE IF THE VECTOR CAN BE RETURNED TO IN M 
STEPSo THE PROGRAM WILL GENERATE DATA DECKS FOR A 
SERIES OF TESTS TO RETURN IN M STEPS-. MSTART IS THE 
STARTING INDEX FOR THE NUMBER OF STEPS TO RETURN, AND 
MSTOP IS THE ENDING INDEXo THESE INDICES ARE READ IN 
AS DATAo THE FIRST VECTOR TO BE TESTED IS CALL FD 
"TEST"o EACH SUCCEEDING VECTOR IS CALLED "REVI" WHERE 
I IS THE NUMBER OF THE REVISION OF THE RIGHT HAND SIDEo 
THE VARIABLE IR IS THE NUMBER OF REVISIONS IN THE 
CONTROL PROGRAMo THE ARRAYS SHOULD BE DIMENSIONED IN 
THE FOLLOWING WAY - 

AX( ( MSTGP-1 »*(N+1 ) + N, ( MSTOP-1 )* (N+l ) ) 

KEY - SAME AS AX 
X ( N » IR + 1 ) 

P(N,N) 

PT ( N , N ) 

PX ( N , N ) 

UNI ( N, N)o 

THE PROGRAM WILL ACCOMODATE A PROMOTION MATRIX WITH 
N = 99 AND WITH MSTOP < 100o 






THERE IS A CERTAIN AMOUNT OF REPETITION BETWEEN THIS 
PROGRAM AND THE FOLLOWING PROGRAM FOR A TEST FOR A 
REACHABLE SET IN M STEPSo TO MINIMIZE DOCUMENTATION, 
CERTAIN MODULES IN THIS PROGRAM NEED ONLY BE REPLACED 
BY OTHER MODULES TO OBTAIN THE NEW PROGRAMo THESE 
MODULES ARE MARKED WITH A NUMBER IN THE LEFT HAND 
COLUMN WHICH MEANS THEY ARE TO BE REPLACED BY THE 
STATEMENTS WITH THE CORRESPONDING NUMBER IN THE NEXT 
SECTION TO OBTAIN THE NEW PROGR AM 0 IN ADDITION, THERE 
ARE SEVERAL STATEMENTS WHICH NEED ONLY BE REMOVED, 
THESE ARE MARKED WITH AN R IN THE LEFT HAND COLUMN, 



INTEGER *2 KEY(111,98) 

DIMENSION AX(111,98), X(13,4), P(13,13), PT(13,13I, 
1UNI( 13,13) , PX(13, 13) 

READ INPUT DATA 

1 READ ( 5 , 600 ) MSTART, MSTOP, N, IR 

600 FORMAT (415 ) 
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CL CL CMCVJ <\J CV] CNJ (SI C\J (V C\J CL CL CL CL CL CL CL CL CL CO CO CO CO fO fO CO 



IRl = IR + 1 

READ (5*601) ( ( P< I , J ) , J = 1 ,N ) , 1=1, N) 

601 FORMAT ( 8 FI Oo 0 ) 

READ (5, 601 ) ( (X( I, J) , I=1,N1 , J=l, IRl) 

MBEG = MSTART - 1 
MEND = MSTOP - 1 
DO 999 M = MBEG, MEND 

INITIAL IZE 

ITEST = 1 
JLIM = M*N + M 
ILIM = JLI K + N 
DO 10 J = 1 * JLIM 
DO 10 I =1 , ILIM 
AX ( I ,J) = Oo 
10 KEY ( I , J ) = C 

CALCULATE THE TRANSPOSE OF THE PROMOTION MATRIX, PT , 

AND LOAD THE IDENTITY MATRIX UNI 

DO 20 1=1, N 
DO 20 J= 1 , N 
UNI ( I, J ) = Oo 

pt ( i , j ) = p< j, n 

IFUoEOoJJ UNI ( I , J ) = lo 
20 CONTINUE 

CALCULATE THE RHS FOR THE FIRST N ROWS 

DO 1 J = 1 » I R1 
DO 2 I = 1,N 
SUM = Oo 
DO 3 K = 1,N 

3 SUM = SUM + PT ( I ,K ) * X(K,J) 

2 PX( I ,J ) = SUM 
1 CONTINUE 

INSERT UNI IN THE FIRST N ROWS AND COLUMNS OF THE 
CONSTRAINT MATRIX, AX 

DO 30 1=1, N 
DO 30 J = 1 , N 
AX ( I , J I = UNI( I , JJ 
30 IF(AX( I , JIoNEoOl KEY ( I , J ) = 1 

INSERT PT AND -I IN THE APPROPRIATE COLUMNS OF THE NEXT 
4 (M-2)*N ROWS OF AX 

LASTC = N 
NXTR = N+l 
NXTC = N+l 
JLASTC = 1 
LASTR = 2*N 
JNXTC = 2* N 
MOD = 1 

5 I F ( MODo E Qo M I GO TO 100 

DO 50 I = NXTR, LASTR 

00 40 J = JLASTC, LASTC 

AX ( I , J ) = PT( I-M0D#N,J-(M0D-1)*N) 

40 I F ( AX ( I , JJoNEoO) KE Y ( I , J ) = 1 

DO 50 J = NXTC, JNXTC 

1 F ( I oEOo J I AX ( I , J ) = - loO 

50 I F ( A X ( I ,J)oEQo-loO) KE Y ( I , J ) = 1 

NXTC = NXTC + N 
NXTR = NXTR + N 
LASTC = JNXTC 
JLASTC = NXTC - N 
LASTR = NXTR + N - l 
JNXTC = LASTC + N 
MOD = MOD + 1 
GO TO 5 
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f 



100 

60 



65 

70 



700 

701 

702 



703 



350 

753 

200 



351 

741 

352 

742 

353 

743 

354 

744 
201 



705 

320 

720 

202 



INSERT PT IN THE CORRECT COLUMNS OF THE NEXT N ROWS OF 

AX 

DO 60 I = NXTR » L ASTR 

DO 60 J = J LASTC , L ASTC 

AX ( I , J ) = PT(I-MOD*N, J-(M0D-1>*N) 

IF(AX( I , JloNEoO) KEY(I»J) = 1 

INSERT l'S IN THE CORRECT COLUMNS OF THF LAST M ROWS OF 

AX WHICH CORRESPOND TO THE FEASIBILITY CONSTRAINT 

MOD = MOD + 1 
IM = LASTR + 1 
IMM = LASTR + M 
DO 70 I = 1 M , I MM 
ICOL = ( I-IM)*N + 1 
JCOL = ICOL - 1 + N 
DO 65 J = ICOL, JCOL 
AX ( I , J ) = lo 

KEY ( I , J ) = 1 
K = M*N + I - LASTR 
AX ( I , K ) = lo 

KE Y ( I , K ) = 1 
CONTINUE 

PUNCH OUT THE DATA SET NAME AND ROWS CONSTRAINTS CARDS 
WRITE(6,700 ) 

FORMAT ( 'NAME* ,T15, ' TEST' ) 

WRITE<6,701 ) 

FORMAT ( • ROWS’ ) 

WRITF( 6,702) 

FORMAT (T2, • N OBJF « ) 

DO 200 J=1 , N 
IF(JoGT 0 9) GO TO 350 
WRITE! 6 ,703 ) J 
FORMAT ( T2 , ' G AO, 'll) 

GO TO 200 
WRITE(6,753) J 
FORMAT ( T2 , ' G A0,',I2) 

CONTINUE 
DO 201 I = 1 , M 
DO 201 J = 1 , N 

IF! IoLEo9o ANDo JoLEoP) GO TO 351 
IF(IoLEo9oANDo JoGT 0 9) GO To 352 
IF!IoGTo9oANDo JoLEo9) GO TO 353 
IF! IoGTo9o AND® JoGT 0 9 ) GO TO 354 
WRI T E ( 6 , 741 ) I,J 
FORMAT ( T2, ' L A', II, ',',11) 

GO TO 201 

WRITE ( 6 , 742 ) I,J 

FORMAT ( T2, ' L A', II, ',',12) 

GO TO 201 

WRI TE ( 6 , 743 ) I,J 

FORMAT ( T 2, ' L A', 12, ',',11) 

GO TO 201 

WRITE! 6, 744 ) I,J 

FORMAT ( T2, ’ L A', 12, ',',12) 

CONTINUE 

DO 202 I = 1,M 

I F ( I o GT o 9 ) GO TO 320 

WRITE! 6 , 705 ) I 

FORM AT ( T2, 1 E AR ' , 1 1 ) 

GO TO 202 

WR I TE ( 6 , 72 0 ) I 

FORMAT ( T2, ' E AR',12) 

CONTINUE 

PUNCH OUT THE COLUMNS DATA CARDS . 

WR ITE ( 6 , 706 ) 
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706 FORMAT ( 'COLUMNS' ) 

JX = 0 
JEND = M*N 
DO 300 J = 1 1 J END 
JX = JX + 1 
DO 302 K = 1,M 

I F ( J 0 LEoK*N, ANDo JoGTo (K-1)*N) IX = K-l 
IF( JoEQo IX*N+1 ) JX = 1 

302 CONTINUE 

DO 301 I = 1, I L I M 
I F ( KEY ( I » J ) o NEo 1 ) GO TO 301 
M2 = M + 2 
DO 303 L = 1 » M2 

IF( IoLEoL^NoANDo IoGTo (L-1)*N) IA = L-l 

303 CONTINUE 

JA = I - IA*N 

IF(IoGTo(M+l)*N> JA = I - ( M+l )#N 

IF(lXoLEo9oANDoJXoLEo9oANDoIA 0 LEoQ,ANDo JAoLEo9) GO TO 

1360 

IF(IXoLEo9oANDoJXoGTo9oANDoIAoLEo9,ANDoJAoLEo9) GO TO 

1361 

IF(IXoLEo9oANDoJXoGTo9oANDoIAoLEo9oANDoJAoGTo9) GO TO 

1362 

IF( IXoLEo9oANDoJXoGTo9oANDoIAoGTo9oANDo JAo LEo9) GO TO 

1363 

IF( IXoLEo9oANDoJX 0 GTo9oANDo IA 0 GTo9oANDo JAoGTo9) GO TO 

1364 

IF( IXoGTo9 0 ANDoJXoGTo9oANDoIAoLEo9oANDo JAoGTo9) GO TO 

1365 

IF( IXoGTo9oANDoJXoGToQoANDoIAoGTo°oANDo JAo LEo9I GO TO 

1366 

IF( IXoGT 0 9oANDoJXoGTo9oANDoIAoGTo9oANDo JAoGTo9) GO TO 

1367 

IF( IXoLEo9o ANDoJXo LEo9oANDoI Ao LEo9oANDo JAo GTo9) GO TO 

1368 

IF( IXoGT 0 9o ANDo JX 0 LEo9o ANDo I AoGTo9o ANDo JAo LE, 9) GO TO 

1369 

IF(IXoLEo9oANDo JX 0 LEo9o ANDo IAoGTo9oANDo JAo GTo 9 ) GO TO 

1370 

IF(IXoGTo9oANDoJXoLEo9oANDoIAoGTo9oANDoJAoGTo9) GO TO 

1371 

IF( IXoLEo9oANDoJXoLEo9oANDoIAoGTo9o ANDo JAo LEo9) GO TO 

1372 

1F(IXoGTo9oANDo JX oLEo9oANDoIAo LEo9oANDo JAoGTo9) GO TO 

1373 

IF(IXoGTo9oANDoJXoGT 0 9oANDoIAoLEo9oANDoJAoLEo9) GO TO 

1374 

IF(IXoGTo9oANDoJXoLEo9oANDoIAoLEo9oANDoJAoLEo9) GO TO 

1375 

360 IFUoGTo (M + l )*N) GO TO 380 
WRITE (6,760) I X, JX , I A, JA , AX! I , J ) 

760 FORMAT ( T5, • X* , 1 1 , • , • , 1 1 , T1 5 , • A • , 1 1 , • , • , 1 1 , T25 , FI Oo ~f ) 
GO TO 301 

380 WRITE (6 , 780 ) IX, JX , JA, AX( I , J ) 

780 FORMAT (T5, ' X • , 1 1 , • , ' , I 1 , T1 5 , • AR' , 1 1 , T25 , FlOo 7 ) 

GO TO 301 

361 IF( IoGTo (M+l )*N) GO TO 381 
WRITE! 6, 761 ) IX, JX, IA, JA, AX( I , J) 

761 FORMAT (T5, • X • , 1 1 , » , • , I 2 , T1 5 , • A • , 1 1 , • , • , I 1 , T2 5 , FI Oo 7 ) 
GO TO 301 

381 WRITE(6,781) IX, JX , JA , AX ( I , J ) 

781 FORMAT ( T5 , • X • , 1 1 , • , • , I 2 , T1 5 , • AR • , 1 1 , T2 5 , F 1 Oo 7 ) 

GO TO 301 

362 I F ( IoGTo(M + l )*N) GO TO 382 
WRITE(6,762) IX, JX , I A, JA , AX ( I , J ) 

762 FORMAT! T5, • X' , 1 1 , • , • , I 2 , T1 5 , • A • , 1 1 , • , • , 1 2 , T2 5 , FlOo 7 ) 
GO TO 301 

382 WRITE(6,782) IX, JX, JA, AX ( I , J ) 

78 2 FORMAT ( T5, • X' , 1 1 , • , • , I 2 , T1 5 , • AR * , I 2 , T2 5 , FI Oo 7 ) 

GO TO 301 

363 IF(IoGTo(M+l)*N) GO TO 383 
WRITE! 6,763) IX, JX , I A, JA , AX! I , J ) 
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763 

383 

783 

364 

764 

3e4 

784 

365 

765 

385 

785 

366 

766 

386 

786 

367 

767 

387 

787 

368 

768 

388 

788 

369 

769 

389 

789 

370 

770 

390 

790 

371 

771 

391 

791 

372 

772 

39 2 

792 

373 

773 



FORMAT ( T5, • X' 


,11 


, ' , 


' , I 2, T1 5 , 


•A* 


,12, 


' , * , 


11, 


T2 5, 


FlOo 


7) 


GO TO 301 






















WRITE(6,783) 


I x, 


JX, 


JA , AX ( I , J 
















FORMAT ( T5 * 'X* 


, 11 


, ' , 


• , I 2, T 1 5, 


• AR 


',11 


, T25 


, F 1 Oo 7 ) 






GO TO 301 






















I F ( IoGTo (M+l ) 


*N) 


GO 


TO 384 
















WR I T E ( 6 » 76 4 ) 


IX, 


JX, 


IA, JA, AX( 


I , J 














FORMAT ( T 5 , ' X* 


,11 


,• , 


' , I 2 , T1 5 , 


•A' 


,12, 


' , • , 


12, 


T25, 


FlOo 


7 ) 


GO TO 301 






















WR ITE ( 6 , 784 ) 


IX, 


JX, 


JA, AX ( I , J 
















FORMAT ( T5 * • X' 


,11 


,' , 


• , I 2 , T1 5 , 


• AR 


' ,12 


, T25 


, FI Oo 7 ) 






GO TO 301 






















I F ( IoGTo (M+l ) 


*N) 


GO 


TO 385 
















WRITE( 6, 765 ) 


IX, 


JX, 


I A, JA ,AX( 


I , J 














FORMAT ( T 5 » »X’ 


,12 


,' , 


' , I 2 , T1 5 , 


•A’ 


, 11 , 


’ , ' , 


12, 


T25, 


FlOo 


7) 


GO TO 301 






















WRITE( 6, 785 ) 


IX, 


JX, 


JA, AX( I ,J 
















FORMAT (T5, • X' 


, 12 


, • , 


' , I 2, T1 5, 


* AR 


' ,12 


, T2 5 


, FI Co 7) 






GO TO 301 






















I F ( IoGTo (M + l ) 


*N) 


GO 


TO 386 
















WRITE( 6, 766 ) 


IX, 


JX , 


IA, JA, AX( 


I , J 














FORMAT ( T5 , • X* 


,12 


,’ , 


• , I 2 , T1 5 , 


•A’ 


,12, 


’ , ' , 


11, 


T25, 


FlOo 


7) 


GO TO 301 






















WRITE(6,786) 


IX, 


JX, 


JA, AX( I , J 
















FORMAT ( T5 , • X* 


,12 


, • , 


• ,1 2,T1 5, 


• AR 


• ,11 


,T25 


, FlOo 7 ) 






GO TO 301 






















I F( IoGTo (M+l ) 


*N) 


GO 


TO 387 
















WRITE( 6,767) 


IX, 


JX, 


IA, JA, AX( 


I , J 














FORMAT ( T5, • X* 


,12 


, • , 


' , I2,T15, 


• A • 


, 12 , 


' , ’ , 


12, 


T25, 


FlOo 


7) 


GO TO 301 






















WRITE (6,787) 


IX 


, JX 


, JA , AX ( I , 


J) 














FORMAT ( T5 , • X» 


, 12 


,' , 


' , I2,T15, 


•AR 


• ,12 


,T25 


, Fl Co 7) 






GO TO 301 






















I F ( IoGTo (M+l ) 


*N) 


GO 


TO 388 
















WRITE ( 6 , 76 8 ) 


IX, 


JX , 


IA, JA, AX( 


I , J 














FORMAT ( T5 , * X* 


,11 


, 


' , I 1 , T1 5 , 


• A* 


,11 , 


• , ’ , 


12, 


T25, 


FlOo 


7 ) 


GO TO 301 






















WR IT E ( 6 , 78 8 ) 


IX, 


JX, 


JA, AX ( I , J 
















FORMAT ( T5 , 'X' 


,11 


, ' , 


' ,I1,T15, 


• AR 


’ ,12 


, T2 5 


, FI Oo 7 ) 






GO TO 301 






















IFUoGTo (M + l ) 


*N) 


GO 


TO 389 
















WRIT E ( 6,769) 


IX, 


JX, 


IA, JA, AX( 


I , J 














FORMAT (T5, * X* 


, 12 


,’ , 


' , I1,T15, 


• A • 


t 12, 


', •, 


lit 


T25 , 


FlOo 


7 ) 


GO TO 301 






















WRITE( 6,789 ) 


IX, 


JX, 


JA, AX ( I ,J 
















FORMAT (T5, ’X* 


,12 


, ' , 


• , I1,T15, 


• AR 


',11 


, T25 


, FlOo 7) 






GO TO 301 






















I F ( IoGTo (M + l ) 


*N) 


GO 


TO 390 
















WRITE(6,770) 


IX, 


JX, 


IA, JA, AX( 


I , J 














FORMAT ( T5, • X' 


,11 


,' , 


' , I 1 , T1 5 , 


•A* 


,12 , 


' , • , 


12, 


T25, 


FlOo 


7) 


GO TO 301 






















WRITE( 6,790 ) 


IX, 


JX, 


JA, AX ( I , J 
















FORM AT ( T5 , 'X* 


, 11 


,’ , 


•,11,115, 


• AR 


• ,12 


, T2 5 


, FlOo 7 ) 






GO TO 301 






















I F ( IoGTo (M+l ) 


*N) 


GO 


TO 391 
















WRITE (6,771 ) 


IX, 


JX, 


IA, JA,AX( 


I , J 














FORMAT (T5, * X* 


,12 


,* , 


• , I1,T15 , 


• A • 


,12, 


' , ’ , 


12, 


T25, 


FlOo 


7 ) 


GO TO 301 






















WRITE(6,791 ) 


IX, 


JX, 


JA, AX( I , J 
















FORMAT (T5, ' X» 


,12 


,• , 


' , I1,T15, 


• AR 


• ,12 


,T2 5 


, FlOo 7) 






GO TO 301 






















I F ( IoGTo (M+l ) 


*N) 


GO 


TO 392 
















WRITE( 6,772 ) 


IX, 


JX, 


IA, JA,AX( 


I , J 














FORMAT (T5, • X' 


,11 


, • , 


• ,I1,T15, 


' A • 


,12 , 


1 t * t 


11 , 


T25, 


FlOo 


7) 


GO TO 301 






















WRITE( 6,792 ) 
FORM AT ( T 5 , • X* 


III 


JX , 
, • , 


JA, AX( I , J 
' , I 1 , T1 5 , 


' AR 


• ,11 


,T25 


, F 1 Oo 7 ) 






GO TO 301 






















I F ( I oGTo ( M+ 1 ) 


*N) 


GO 


TO 393 
















WRITE( 6,773 ) 


IX, 


JX, 


IA, JA,AX( 


I , J 














FORMAT (T5, ' X* 


,12 




* , 1 1 , T1 5 , 


•A* 


,11 , 


' , • , 


12, 


T25, 


FlOo 


7 ) 


GO TO 301 
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393 WR I T E ( 6 , 793 ) I X , J X , J A , AX ( I , J ) 

793 FORMAT (T5, 1 X' , 12 , • , • , I1,T15, • AR' , I 2,T2 5, FlOo 7) 

GO TO 301 

37 A IF(IoG T o(M + l)*N) GO TO 394 

WR I T E ( 6 , 77 4 ) I X , J X , I A , JA , AX ( I , J ) 

774 FORMAT (T5, 'X',I2,',',I2,T15,'A',Il,',',Il,T25,F10o7) 

GO TO 301 

394 WR I TE ( 6 , 794 ) I X , JX , JA , AX ( I , J ) 

794 FORMAT (T5, • X • , I 2 , • , M 2 , T1 5 , ' AR • , 1 1 , T25 , F 10o 7 ) 

GO TO 301 

375 I F ( IoGTo (M+l )*N) GOTO 395 

WR I T E ( 6 , 77 5 ) I X , J X , I A , J A , AX ( I , J ) 

775 FORMAT (T5, ' X • , I 2 , • , • , 1 1 , T1 5 , ' A • , 1 1 , • , ' , 1 1 , T25 ♦ FlOo 7 ) 

GO TO 301 

395 WR I T E ( 6 * 79 5 ) I X, JX , JA , AX( I , J ) 

79 5 FORMAT (T5, • X • , 1 2 , • , ' , 1 1 , T1 5 , ' AR • , 1 1 , T 25 , F 1 Oo 7 ) 

301 CONTINUE 

I F < JoLTo JEND) GO TO 300 
DO 305 IZ = 1 » M 
JZ = IZ 

I F ( I Zo GT 0 9 ) GO TO 304 
WRITE! 6 , 709 ) IZ, JZ 

709 FORMAT (T5, • XA • , 1 1 , T15 , • OBJ F • , T2 5 , • lo ' ,T40, 'AR* ,11,^50, 

l'lo' ) 

GO TO 305 

304 WR I T E ( 6 , 70 8 ) IZ, JZ 

708 FORMAT (T5, • X A • , I 2, T 15 , • OBJ F • ,T 25, • lo ' , T40 , ' AR • , I 2 , T50 , 

l'lo') 

305 CONTINUE 
300 CONTINUE 

PUNCH OUT THE RIGHT HAND SIDE DATA CARDS 

WRITE! 6,710) 

710 FORMAT! 'RHS' ) 

DO 325 I = 1 , N 

I F ( I oGT 0 9 ) GO TO 333 
WRITE! 6,717) I, PX ( I , 1 ) 

717 FORMAT (T5, • B1 • , T15 , • AO, ' ♦ 1 1 ,T2 5, F10 o 7) 

GO TO 325 

333 WRITE! 6,718) I, PX ( I , 1 ) 

718 FORMAT <T5, • El • , T1 5 , • AO , ' , I 2 , T2 5, FI Oo 7 ) 

325 CONTINUE 

DO 306 I » 1, N 

IF(MoLEo9oANDoIoLEo9) GO TO 311 
IF(MoLEo9oANDoIoGTo9) GO TO 312 
IF(MoGTo9o ANDo IoLEo9) GO TO 313 
IF(MoGTo9oANDoIoGT 0 9) GO TO 314 

311 WR I T E( 6 , 71 1 ) M, I, X ( I ,1) 

711 FORMAT (T5, • B1 ' , T1 5 , • A » , 1 1, • , • , 1 1 ,T 25 , FlOo 7 ) 

GO TO 306 

312 WR IT E ( 6 , 72 2 ) M, I, X(I,1) 

722 FORMAT (T5, 'B1',T15,'A»,I1,',', I 2 ,T25 , FlOo 7 ) 

GO TO 306 

313 WRITE(6,723) M, I, X ( I , 1 ) 

72 3 FORMAT (T 5, • B1 ' , T 1 5 , • A • , I 2 , ' , • , 1 1 ,T25 , FlOo 7 ) 

GO TO 306 

314 WRITE! 6 , 724 ) M, I, X ( I , 1 ) 

7 24 FORMAT !T5, ' B1 ' »T15 , ' A ' ,12,',', 1 2 ,T2 5 , FlOo 7 ) 

306 CONTINUE 

DO 307 I = 1 * M 
I F ( I oGTo Q ) GO TO 315 
WRI T E ( 6 , 71 2 ) I 

712 FORMAT (T5, • El • ,T1 5 , • AR • , 1 1 , T25 , ' lo ' ) 

GO TO 307 

315 WR ITE! 6 , 72 5 ) I 

72 5 FORMAT (T 5, ' B 1 ' , T 1 5 , • AR • , 1 2 , T25 , ' lo ' ) 

307 CONTINUE 

WRI T E ( 6 , 71 3 ) 

713 FORMAT! • ENCATA' ) 

801 IF! ITESToEQo IR1 ) GO TO 999 
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PUNCH OUT THE APPROPRIATE RHS REVISION DATA CARDS 

I F ( I TESTo GT 0 9 ) GO TO 308 
WRITE(6,714) ITEST 

714 FORM AT ('NAME* , T1 5, ' REV ' , 1 1 ) 

GO TO 30R 

308 WRITE(6,715) ITEST 

715 FORMAT! 'NAME' ,T15, 'REV* ,12) 

309 WRITE(6,71C) 

WRITE! 6, 71 6) 

716 F0RMAT(T3, 'MODIFY* ) 

I T 1 = ITEST + 1 
DO 326 I = 1 , N 
IF(IoGT 0 9) GO TO 334 
WR I T E ! 6 , 71 7 ) I,PX(I,IT1) 

GO TO 326 

334 WR I T E( 6 , 71 8 ) I, PX(I,IT1) 

326 CONTINUE 

DO 310 I = 1, N 

IF(MoLEo9o ANDo IoLEo9) GO TO 321 
IF!MoLEo9o ANDo IoGT 0 9) GO TO 322 
IF!MoGTo9oANDo IoLE 0 9) GO TO 323 
IF(MoGTo9o ANDo IoGTo9) GO to 324 

321 WR I T E ( 6 , 71 1 ) M, I, X(I,IT1) 

GO TO 310 

322 WRITE<6,722) M, I, X(I,IT1) 

GO TO 310 

323 WR IT E ( 6 , 72 3 ) M, I, X(I,IT1) 

GO TO 310 

324 WR IT E ( 6 , 724 ) M, I, X(I,IT1) 

310 CONTINUE 

WR ITE < 6 , 71 3 ) 

ITEST = ITEST + 1 
GO TO 801 
999 CONTINUE 
STOP 
END 
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****************************** 



* * DATA GENERATOR - TEST FOR A REACHABLE SET IN M STEPS * * 



****** ROBERT Lo ARMACOST - - - 1 JULY 1970 ****** 



************* ***************** 



* * * * _* ************************* 

THIS PROGRAM WILL PUNCH OUT A SET OF DATA CARDS FOR A 
LINEAR PROGRAM USING THE MPS PACKAGE ON THE IBM 360/67, 
THE BASIC INPUT TO THE PROGRAM IS A SQUARE PROMOTION 
MATRIX OF SIZE N, AND K VECTORS WHICH ARE TO BE TESTED 
TO DETERMINE IF THEY ARE ELEMENTS OF THE SET B<M>, 

THE PROGRAM WILL GENERATE DATA DECKS FOR A SERIES OF 
REACHABLE SETSo MB EG IS THE STARTING INDEX FOR THE 
SUBSCRIPT OF THE TESTo MEND IS THE ENDING INDEX FOR 
THE SUBSCRIPT OF THE TESTo THESE INDICES ARE READ IN 
AS DATAo THE FIRST VECTOR TO BE TESTED IS CALLED 
"TEST" 0 EACH SUCCEEDING VECTOR IS CALLED "REVI" WHERE 
I IS THE NUMBER OF THE REVISION OF THE RIGHT HAND SIDEo 
THE VARIABLE IR IS THE NUMBER OF REVISIONS IN THE 
CONTROL PROGRAMo THE ARRAYS SHOULD BE DIMENSIONED IN 
THE FOLLOWING WAY - 

AX( (MENC)^(N + 1 ) + N, MEND* ( N + l ) ) 

KEY - SAME AS AX 
X(N, IR+1 ) 

P(N,N) 

PT(N,N) 

PTKN,N)« 

THE PROGRAM WILL ACCOMODATE A PROMOTION MATRIX WITH 
N = 99 AND WITH MEND < 100o 



****************************** 



****************************** 

THIS PROGRAM IS CONSTRUCTED BY REMOVING THOSE CARDS 
MARKED WITH AN R IN PROGRAM TO TEST FOR RETURN IN M 
STEPS, AND TC REPLACE THOSE MODULES MARKED WITH NUMBERS 
BY THE MODULES LISTED BELOW WITH THE CORRESPONDING 
NUMBERSo 

****************************** 



1 RE AD ( 5 , 600 ) MB EG, MEND, N, IR 

CALCULATE THE TRANSPOSE OF THE PROMOTION MATRIX, PT, 
AND CALCULATE I - PT WHERE I IS THE IDENTITY MATRIX 

DO 20 1=1, N 
DO 20 J = 1 , N 
PT( I ,J » = P < J , I ) 

PTI ( I, J) = -PT ( I , J ) 

IF(IoEOoJ) PT I ( I , J ) = lo + PT I ( I , J ) 

20 CONTINUE 

INSERT PTI IN THE FIRST N ROWS AND COLUMNS OF THE 
CONSTRAINT MATRIX, AX 
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DO 30 1=1, N 
DO 30 J = 1 » N 
AX< I ,J ) = PTI ( I , J) 

30 IF( AX( I , J)o NEoO) K EY ( I » J ) = 1 

4 ( M- 1 » ROWS OF AX 
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****************************** 



* STEADY STATE SET - EXTREME POINTS AND TEST FOR INCLUSION * 



****** ROBERT Lo ARMACOST - - - 1 JULY 1970 ****** 
****************************** 



THE MAIN PROGRAM IS DIMENSIONED TO HANDLE A PROMOTION 
MATRIX WITH 15 STATES, THE PROMOTION MATRIX IS P AND 
ITS INVERSE WHICH HAS BEEN NORMALIZED IS PIN- THE 
ROWS OF PIN ARE THE EXTREME POINTS OF THE STEADY STATE 
SETo 

THE EXTREME POINTS OF THE STEADY STATE SET ARE USED 
TO DETERMINE IF A TEST POINT IS IN THE STEADY STATE 
SET BY CALCULATING THE NECESSARY COEFFICIENTS TO 
MAKE THE TEST POINT A LINEAR COMBINATION OF THE 
EXTREME POINTSo IF THE SUM OF THE COEFFICIENTS IS 
GREATER THAN ONE, OR IF ANY COEFFICIENT IS NEGATIVE, 

WE CONCLUDED THAT THE TEST POINT IS NOT IN THE STEADY 
STATE SETo THE COEFFICIENTS ARE DENOTED ALF, AND THE 
TEST POINTS APE X., THE PROGRAM IS SET UP TO ACCEPT, 
FIVE SUCH TEST VECTORS WITH FIFTEEN RANKS- 

THE SIZE OF THE PROMOTION MATRIX AND THE LENGTH OF THE 
TEST VECTOR ARE INPUT TO THE PROGRAM IN THE VARIABLE 
NNo THE NUMBER OF TEST VECTORS IS INPUT IN THE 
VARIABLE KKo 






DIMENSION UNI ( 15 , 1 5 ) » P ( 1 5 , 1 5 ) , PIN<15,15), L(15), 

1M( 15 ) , SUM (151, ALF ( 15 I * X ( 1 5 , 5 1 
RE AD ( 5 » 599 ) NN, KK 
599 FORMAT (215) 

CALL WORM UNI , P , PI N, L , M , AL F , X , NN , KK , SUM ) 

CALL EXIT 
END 



****************************** 

* * * SUBROUTINE WORK - THIS SUBROUTINE CALCULATES THE * * * 

* * * EXTREME POINTS OF THE STEADY STATE SET, M, AND * * * 

* * * TESTS EACH TEST VECTOR FOR INCLUSION IN THAT SET * * * 



SUBROUTINE WORK ( UN I , P , P I N , L , M , ALF , X , NN , KK, SUM) 
DIMENSION UNI (MN,NN ) , P ( NN, NN ) , PIN( NN,NN) , L (NN) , 
1M(NN), SUM ( NN ) , AL F ( NN ) , X ( NN , KK ) 

RE AD (5, 600 ) <(P(I,J),J=1,NN),I=1,NN) 

6C0 FORMAT ( 8F1 Oo 0 ) 

RE AD (5, 601 ) <<X(J,K),J=1,NN),K=1,KK) 

601 FORM AT ( 8 FIO 0 O) 

INITIALIZE 

DO 10 I = 1 , NN 
L ( I ) = 0 
M( I ) = 0 
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SUM! I) = Oo 

DO 10 J = 1,NN 

UNI ( I, J ) = - P ! I , J ) 

IF(IoEOoJ) UN I < I , J ) = lo + UNI(I,J) 

10 CONTINUE 

WRITE! 6,699 ) 

WRITE OUT PROMOTION MATRIX 
DO 2 1=1, NN 

WR I T E ( 6 , 700 ) ( P( I , J I , J=1 ,NN) 

WR IT E ( 6 , 71 0 ) 

2 CONTINUE 

CALCULATE THE INVERSE OF P 

CALL MINV(UM,NN,D,L,M) 

CALCULATE THE NORMALIZED INVERSE OF P, PIN 

DO 20 1=1, NN 
DO 15 J = 1 » NN 

SUM! I ) = SUM! I ) + UNI ( I , J) 

15 CONTINUE 

DO 25 J=1 , NN 

PIN! I , J ) = UNI ( I , J ) / SUM! I ) 

25 CONTINUE 
20 CONTINUE 

WRITE! 6,699) 

699 FORMAT ( 1H1 ) 

WRITE OUT NORMALIZED INVERSE OF P 
DO 5 1=1, NN 

WRITE! 6, 700 ) ( PI N( I , J ) , J=1 ,NN) 

WRITE! 6,710) 

5 CONTINUE 

700 FORMAT ( 1 5F 8o 6 ) 

710 FORM AT ( 1 H ) 

TEST TO SEE IF TEST VECTOR IS IN THE STEADY STATE SET 

DO 130 K=1 , KK 
MM=K 

DO 1 1=1, NN 
1 ALF! I ) = Oo 

ALF(l) = X ( 1 , K ) / PI N ( 1 , 1 ) 

TEST = ALF ( 1 ) 

JEND = NN- 1 
DO 120 J =2 , JEND 
I END = J - 1 
SUMD = Oo 
DO 110 I =1 , I END 

110 SUMD = SUMD + ALF! I ) * P I N ( I , J ) 

ALF! J ) = ( X ( J , K ) - SUMD) / PIN(J,J) 

TEST = TEST + ALF(J) 

IF(ALFU)oLToO) GO TO 100 
IF(TEST 0 GTol ) GO TO 200 
120 CONTINUE 

ALF( NN ) = lo - TEST 
I F( ALF ( NN) o LTo 0) GO TO 100 
IF(ALF(NN)oLToX(NN,K) ) GO TO 200 
130 CONTINUE 

WRI TE ( 6 , 702 ) 

702 FORMAT (1H1 ,T20, ’THE VECTOR IS IN THE STEADY STATE SET' 
1//T2C, ' ALF COEFFICIENTS' ,T45, 'TEST VECTOR'//) 

DO 201 LL= 1 , NN 

WRITE(6,703) ALF(LL), X ( LL , MM ) 

703 FORMAT (T20, F10o8,T45,F10 o 8 ) 

201 CONTINUE 

IF(KoLToKK) GO TO 130 
GO TO 999 
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100 WRITE<6,704) 

704 FORMAT ( 1H1 ,T20 ,' THE VECTOR IS NOT IN THE STEADY STATE' 
It* SET, ALF IS NEGATIVE' //T20, 'ALF COE FF I C I E NTS • , T45 , 

2 ' TEST VECTOR'//) 

DO 202 LL= 1 * NN 

WRITE! 6,703) ALF(LL), X ( LL ,MM ) 

202 CONTINUE 
IF(KoLToKK) GO TO 130 
GO TO 99Q 

200 WRITE(6,705) 

705 FORMAT (IH1 ,T20, ' THE VECTOR IS NOT IN THE STEADY STATE' 
1,’ SET, TEST IS GREATER THAN ONE'//T20, 

2 ' ALF COEFFICIENTS' » T45 , 'TEST VECTOR'//) 

DO 203 LL= 1 , NN 

WRITE(6,703) ALF(LL), X ( LL , MM ) 

203 CONTINUE 
IF(KoLToKK) GO TO 130 

999 RETURN 
END 
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